Antiferromagnetic ordering of energy levels for spin ladder with four-spin cyclic 
exchange: Generalization of the Lieb-Mattis theorem 
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The Lieb-Mattis theorem is generalized to an antiferromagnetic spin-ladder model with four-spin 
cyclic exchange interaction. We prove that for J > 2K, the antiferromagnetic ordering of energy 
levels takes place separately in two sectors, which remain symmetric and antisymmetric under the 
reflection with respect to the longitudinal axis of the ladder. We prove also that at the self-dual 
point J = 2K, the Lieb-Mattis rule holds in the sectors with fixed number of rung singlets. In both 
cases, it agrees with the similar rule for Haldane chain with appropriate spin number. 



The systems with multi-spin exchanges have gained a 
lot of interest for a long time (for a recent review, see 
Refs. [Hill). These interactions arise at higher orders of a 
strong coupling expansion of a half-filled Hubbard model 
and provide perturbative corrections to the Heisenberg 
model.— The multi-spin cyclic permutations first were 
suggested to be important in two-dimensional (2D) quan- 
tum solids like 3 He4i The relevance of the four-spin cyclic 
interaction around square plaquettes of CuC>2 planes in 
high-temperature superconductors was suggested in Ref. 
[H, and then it was proven experimentally^* Recently, their 
quasi- ID counterparts with similar structure and prop- 
erties have been studied intensively (see Ref. for a 
review). Note that the spin-ladder model is the simplest 
system, where the four-spin cyclic exchange appears from 
the electron interaction. In fact, the investigations of 
some copper-based spin-ladder materials have revealed 
the relevance of four-spin exchange termi (see Ref. d for 
a review). The ground-state phase diagrams of frustrated 
spin- 1/2 ladder systems are well investigated.— The inclu- 
sion of four-spin interactions may result in new unconven- 
tional phases. 10 ' 11 ^ 2 ' 13 The different phases are related 
by a duality transformation. 12 

In this paper, we will generalize the well-known Lieb- 
Mattis theorem on ordering of energy levels to the ladder 
model with four-spin cyclic interaction. For finite-size 
Heisenberg models on bipartite lattices, Lieb and Mattis 
proved that the lowest energy E(S) in the sector, where 
the total spin is equal to S, is a monotone increasing 
function of the spin for any S > S gs r^ where S gs is the 
spin of the ground state. This property is known as Lieb- 
Mattis theorem about antiferromagnetic ordering of en- 
ergy levels. The bipartiteness means that the lattice can 
be divided into two sublattices A and B, so that all in- 
teractions within the same sublattice are ferromagnetic 
while the interactions between different sublattices are 
antiferromagnetic. Moreover, the quantum ground state 
of finite-size system is a unique multiplet with total spin 
S gs = \Sa — Sb\, which coincides with the spin of the 
classical ground state, namely, the Neel state. Here, Sa 
and Sb are the highest possible spins on corresponding 
sublattices^ In one dimension, the Lieb-Mattis theorem 
is valid for a more general class of quantum systems. In 
particular, it is true for the Hubbard chain. 15 Recently, 



it has been generalized to SU(n) symmetric chainJ^ A 
ferromagnetic ordering of energy levels has also been for- 
mulated and proven for the Heisenberg chain^ 

According to numerical simulations, a weak frustra- 
tion may preserve the antiferromagnetic ordering of en- 
ergy levels^ and the ground-state spin valued whereas a 
stronger frustration can violate the Lieb-Mattis property. 
For many frustrated systems, the lowest levels E(S) show 
approximately parabolic or linear growth^ So, although 
the Lieb-Mattis theorem is not valid for frustrated spin 
systems in general, it (or its proper extensions) may be 
valid for certain systems. In particular, its generalization 
to reflection-symmetric frustrated spin- 1/2 ladder model 
has been formulated and proven recently^ Here, we ob- 
tain similar results for the frustration caused by four-spin 
ring interaction. 

The Hamiltonian of the system reads: 
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H — E J; (Si,; • Si^+i + S 2 ,( • Sa,i+i) 
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where Si j and S2.1 are the spin- 1/2 operators of the first 
and second chains respectively. The cyclic ring exchange 
pU _|_ pU-i j s a superposition of clockwise and counter 
clockwise permutations of four spins around each plaque- 
tte. We consider the following range of couplings: 



jf > 2Ki 



> 0. 



(2) 



The system possesses SU(2) spin symmetry. It has 
also Z2 symmetry corresponding to the reflection with 
respect to the longitudinal axis. So, the Hamiltonian 
remains invariant on individual sectors with fixed values 
of spin S and reflection a = ±1 quantum numbers. 

We will prove that for the model {I), the antiferro- 
magnetic ordering of energy levels holds independently in 
symmetric (<r = 1) and antisymmetric (tr = —1) sectors 
and conforms to the similar rule for Haldane chain,— i.e. 
spin-1 Heisenberg chain, with N and N —1 spins, respec- 
tively. Namely, the lowest-energy levels E a {S) in sectors 



2 



with spin S and parity a are nondegenerate and mono- 
tone increasing functions of 5* for S > S gs (a). Here, 



S gs {a) 



0, if (7 



-1) 



1, ifa=(-l) 



N-l 



is the ground-state spin value in the sector with parity 
a. The nondegeneracy means that all states on corre- 
sponding level form a unique multiplet. So, the ground 
state in a = (— 1)^ sector is a unique singlet while in 
u = (— l)^ -1 sector it is a unique triplet. 

We will prove also that at the self-dual point jf = 
2Ki, the Lieb-Mattis rule holds in the sectors with fixed 
number of rung singlets N and agrees with the similar 
rule for the Haldane chain with N — Nq spins. 

We begin by introducing the basis of three triplet and 
one singlet states for each rung: 



|1> 
|0) 



r 

T 



i 



(3) 



Below, we will prove that all nonvanishing off-diagonal 
elements of the Hamiltonian (HJ,© become negative in 
the basis 



\mi, . . .,m N ) 

= r_l)[JVo/2]+JV o6 +M odd 



mi) 



\m N ) , 



(4) 



where mi = ±1,0,0 marks one of the four rung states 
©. In the sing factor, iVo is the number of singlets, 
and N q is the number of pairs (0, 0) in the sequence 
mi, . . . , tojvj where is on the left-hand side from 0i2i 
A/odd = J2i m 2/-i is the total z-projection of odd-site 
spins. Note that the basic states above are eigenstates of 
the reflection operator with eigenvalue a — (— l) Wo . 

First, we rewrite the Hamiltonian in a form that is 
more convenient for further purposes. The four-spin in- 
teraction term can be expressed via spin operators as 
follows (see, for instance, Ref. l23l ): 



P, 



P 



□-1 



1,1 + 1 ~f r l.l+l — Si, J • Si,z+i + S2,Z • S 2 ,;+l + Si,; ■ 02,1 

+ Si,j+i • S2,;+i + Si,; • S2J+1 + Sijj+i • 02,1 

+ 4(Sij ■ Si,'+i)(S2,i ■ S 2 j+i) 
+ 4(Si^ • S2,;)(Si ! ; + i • S 2 ,'+l) 

- 4(Si, r S 2 ,,+i)(S 2 ,, • Si, m ) + 1/4. 

Using the relation Si • S2 = P12/2 — 1/4, where the op- 
erator P12 permutes two spin states, one can present the 
expression above in the following form: 



+ P U+l - 2 ( S i,i ' S 2.(+l + Si,;+i • S 2 ,i) 



P l,l+i + p l,i+i p u+i- 



(5) 



Here, P}\ +1 -, -P;7+i anc ^ p /i+i are > respectively, four-spin 
permutations along the plaquette rungs, legs and diag- 
onals. Further, we express the two-spin interactions in 



terms of the symmetrized and antisymmetrized rung spin 
operators 



S1.1 + S 



2.;. 



j (a) 



Si. 



»2J- 



(6) 



The operator S\ describes the total spin of lih rung. 
Using ([5]) , @ and omitting nonessential scalar term, one 
can reduce the Hamiltonian (TTJ) to the following form: 



N-l 



1=1 

N-l 



E A ^ F M+i+ p a + i- p u + i) 



N 



l /a(*>)^2 



1 = 1 



In the above equation, we have introduced the sym- 
metrized and antisymmetrized couplings 



Jf 



■J, 



TO 



K,. 



(8) 



Note that a similar decomposition for the Hamiltonian 
without four-spin exchange was applied in Refs. [2H and 



Dermutations P, 



l,l- 



-1 and P* 



have been used 



|13j. Thej 
in Ref. Hj 

The J 1 - part of the Hamiltonian is just the sum of 
rung spins squares, which is diagonal in the basis ([?]). 

The local terms f are also diagonal (with eigenvalues 
±1), since any triplet (singlet) rung state stays symmetric 
(antisymmetric) under the reflection. 

The J s terms correspond to the so-called composite 
spin model.— They conserve the spins of individual rungs 
because S^ describes the total spin of lib. rungi^i The 
singlets remain frozen at their points, and, therefore, the 
factor (— l)[ Ar °/ 2 ]+ Ar oo in (|4|) remains invariant. All non- 
vanishing off-diagonal matrix elements come from the 
exchanges |6) <g> |5) <-> |±1) <g) and |5) (g> |±1) <-> 
|±1) <8> |5) of two neighboring triplet states, which alter 
the sign of (— l) Modd . Note that they coincide with simi- 
lar matrix elements of the Haldane chain, in (nonpositive) 



M„ 



mi ® ... (X) \m N ) 



basis formed by the states (—1) 
So, the composite spin part of the Hamiltonian is non- 
positive in the basis (fj]). 

The matrix elements produced by the antisymmet- 
ric local terms of Hamiltonian |(7J) have been considered 
already in Ref. l2~ll . In terms of lowering-raising op- 
erators S ,< - a - l± = S^ x ± iS^ v , each such term reads 



/ Q (a)+ Q (a)- 
\°l °l + l 



s, 



(«)- 



-,(a)* Q (a)z 
5 Z °l+l ■ 



In contrary 



to the symmetric case, the antisymmetrized spin opera- 
tors mix triplet and singlet states. Their nonzero matrix 
elements are^ 



(0|S (a)+ |-l) = y/2, 
<6|5 (a)z |0) = 1. 



(l|S (a)+ |0) = -y/2, 



(9) 
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Using the above equations and the definition of basic 
states (U]), it is easy to check that 



.,0 t ,0 



1+1, 



q( a ) + c(°)± 
D l °l + l 



-2(-l) [JVn/2l+JV oB+^odd-K/2]-W ( ; 8 -M^ d 



(10) 



where the unchanged sites are replaced by dots. The 
quantum numbers of bra- and -ket states are mentioned, 
respectively, without and with primes. Indeed, according 
to the definition of iV g, the difference N 00 - 
an even number. Also, 7VX = Nn + 2 and M ( 



in (HOD is 



odd 



Kdd±i 



depending on whether I is even or odd. Therefore, the 
exponent in (|10|) is an even number, and the equation is 
true. The next nontrivial matrix element is also negative. 
Namely, 



■ , 0;, ±l/+i, • 



Q (a)=F o(a)± 



■1) 



-M' 



1 



i±l 



(11) 



because the quantum numbers Nqq and No are the same 
for both states. 

In contrast, the z projections of antisymmetrized spin 
operators preserve the quantum number M 0< xd- They 
produce negative matrix elements too: 



,0;,0;+i, . 



(a)z n(a)z 



l+l 



-1) 



Ar oo-K R + (^o-K)/2 



(-1) 



0i,0,+i, 

even— 1 



> 



and 



<...,6,,0j 



q(a)z q(a)z 

°l ^1+1 



(-1) 



-N'- 



,0i,0 H 
(-1) 1 



(12) 



(13) 



The expressions ([ID]). (fTTT) . (fl2"]). (|Tg]) together with con- 
jugate ones constitute the full set of nontrivial matrix 
elements of Hamiltonian generated by J a terms. 

Finally, consider the off-diagonal terms, which are re- 
sponsible for four-spin cyclic exchange. The operator 
Pfl+i j us f permutes two adjacent rung states. At the 
same time, P^ l+ i is a signed permutation: While permut- 
ing the singlet with a triplet state it produces an addi- 
tional minus sign. Therefore, the difference — Pu + i 
vanishes if the spins of both rungs are the same. If their 
spins differ, this operator just permutes them multiply- 
ing by 2. Then, using the definition (H|) of basic states, 
we obtain 



I Pi 



l,U 



2(-l) 



1,1+1 

id-M' , 



■ ■ ,0j,ij+i, 



(14) 



where t = 0, ±1 is any triplet state. Indeed, in the sign 
factor above, M D dd = A/^dd anc l -^oo — = 1 f° r ^ = 6. 
For t = ±1, iV o5 - NL and |M odd -M^ dd | = 1. Together 
with the conjugate element, this is a sole nonvanishing 
off-diagonal matrix element produced by the four-spin 
cyclic exchange term. 

According to the constraints (2|| imposed on the cou- 
plings, the coefficients Jf, J?, Ki in (|7|) are positive. This 



finishes the proof that the ladder Hamiltonian has no pos- 
itive off-diagonal element in the basis (j4|). 

Due to the spin and reflection symmetries, the Hamil- 
tonian is invariant on each (M, a) subspace, all states of 
which have S z — M and a = ±1 quantum numbers. Any 
two basic states within the same subspace are connected 
at least by two-spin interaction terms of the Hamilto- 
nian, as can be easily verified by induction^ So, we can 
apply the Perron-Frobenius theorem^! to the matrix of 
the Hamiltonian restricted to any (M, a) subspace. As 
a result, the lowest energy state there (usually called a 
relative ground state) is unique and is a positive super- 
position of all basic states: 



in) 



M.a 



^ t u mi...tnj 
|mi,...,mjv)G(M,tr) 



\mi, . . .,m N ) , (15) 



where w„ 



> 0. The uniqueness implies that this 



state must have a certain value of spin Sm,<j, which can 
be obtained by comparing with the similar state of the 
Haldain chain. The last model corresponds to the re- 
striction of the composite spin model s[ s ' • S[l\ on 
the states with triplets on all rungs. In fact, all such 
type states from (f4]) form a nonpositive basis for the Hal- 
dane chain.— Its relative ground state |fi)«f in S z = M 
subspace is a positive superposition of all basic states 
and has the highest possible spin value, i.e., \M\, except 
M = and odd N case when its spin is oneJ^ Therefore, 
for (7 = 1, both states \Q)m<j arK l \®)m overlap, and, 
hence, must have the same spin. Similarly, the restric- 
tion of composite spin model to the subspace of states 
with one singlet frozen at the last rung corresponds to 
the Haldane chain with N — 1 sites. So, for a = — 1, the 
spin of (|15p coincides with the spin of the corresponding 
state of the Haldane chain having one site less. Therefore, 
the spin of the relative ground state (fTS")) is \M\ except 
M = and a = (— l)^^ 1 case when it equals one. 

Now we are ready to finish the proof of our main re- 
sult. For a = (— 1) N , the relative ground states |^) ±So . 
are, correspondingly, the highest and lowest states of a 
unique spin-5 multiplet, which has the minimum energy 
E a (S) among all spin-S" states with parity a. The (S, a) 
subspace contains a representative from any multiplet 
with S' > S and parity er. Together with the unique- 
ness condition, this implies that E a (S') must be higher 
than E<j(S) for all S' > S. Consequently, E a (S) is a 
monotone increasing function of S, and the ground state 
in a = (— l) w sector is a nondegenerate spin singlet. For 
a = (— l)^ 1 , the states |fi)±g c are the highest and low- 
est states only if S > 1. Therefore, E a (S) monotonically 
increases in this region. Note that the states |fi) , j and 

|n)oo-> wn i cn have the lowest energy in a = (— l)^ -1 
sector, form a spin triplet. 



Consider separately the limiting case of J" 



2Ki 



when the J a terms in (J7]) are absent. Then 53, (S T ) z 
commutes with the Hamiltonian and with the total spin 
operator^ As a result, the symmetry group SU(2)xi?2 
expands to SU(2)xU(l)=U(2). The U(l) symmetry re- 
flects the invariance under the duality transformation 



4 



and results in the conservation of singlet number NqJ^ 
Any (M, a) subspace splits into invariant subspaces with 
fixed singlet number obeying (— 1)-"° = a. ft is easy 
to see that the Hamiltonian is connected on every such 
subspace. Therefore, the relative ground state there is 
unique and is a positive superposition of all basic states 
(|2]) with S z = M and Nq rung singlets. Again, compar- 
ing it with the action of the composite spin model on 
the states with singlets on the last Nq rungs, one can 
conclude that the relative ground state has the highest 
possible spin value, except M = and odd N — N Q case 
when it is a triplet state. The antiferromagnetic ordering 
of energy levels takes place independently in any sector 
with fixed singlet number and corresponds to the similar 
rule for the Haldane chain with N — Nq spins. 

For appropriate values of couplings, the results of this 



paper remain true, if diagonal interactions Si : ; • S 2i z_|_i + 
• S 2 ,; are included in the Hamiltonian ([T]). Simi- 
larly, one can consider a more general biquadratic inter- 
action of type K((Si,t ■ Si,j +1 )(S a ,i ■ S 2 ,j+i) + 4if i (S 1 , i • 
S2,i)(Si,i+i • S 2 ,(+i) - 4JT{(Si,; • S 2 ,z+i)(S 2 ,; • Su+i) with 
arbitrary couplings K[ because the first term in the sum 
is diagonal in the basis |(3J|, Q. 
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